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ABSTRACT 

We present an improved analytic calculation for the tidal radius of satellites and test 
our results against N-body simulations. 

The tidal radius in general depends upon four factors: the potential of the host 
galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star 
within the satellite. We demonstrate that this last point is critical and suggest using 
three tidal radii to cover the range of orbits of stars within the satellite. In this way 
we show explicitly that prograde star orbits will be more easily stripped than radial 
orbits; while radial orbits are more easily stripped than retrograde ones. This result has 
previously been established by several authors numerically, but can now be understood 
analytically. For point mass, power-law (which includes the isothermal sphere), and a 
restricted class of split power law potentials our solution is fully analytic. For more 
general potentials, we provide an equation which may be rapidly solved numerically. 

Over short times ( ;$ 1 — 2 Gyrs ~ 1 satellite orbit), we find excellent agreement 
between our analytic and numerical models. Over longer times, star orbits within the 
satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this 
causes a convergence of the three limiting tidal radii towards the prograde stripping ra- 
dius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially 
anisotropic. 
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1 INTRODUCTION 

The tidal radius is the radius at which a star within a satel- 
lite becomes unbound or stripped and becomes bound in- 
stead to the host galaxy about which the satellite orbits. The 
problem of calculating the ti dal radius of satellite s has a long 
history which dates back to Ivon Hoernerl (119571) who anal- 
ysed the problem in the context of globular clusters around 
the Milky Way. 

It is a problem with a very wide scope of application, 
from understanding the sharp edge observed in some nearby 
globular clusters iKind ll962'l. to semi-analytic modelling 
of galaxy formation, where interactions and mergers play 

a vital role in our current hierarchical formation paradigm 
j _ II i 

(Taffo ni et al.ll2005fl . It is a problem which has also experi- 
enced a renaissance with the advent of so-called 'near-field 
cosmology'. Our current cosmological models are beginning 
to give precise predictions for the dist ribution, abundanc e 
and morphology of satellite galaxies (iMoore et al.l [1999). 
while data from the nearest satellites within the Local Group 
can p rovide excellent constraints on these theories (iMateol 
1998). Understanding in detail the effects of a tidal field on 



Email: jir22@ast.cam.ac.uk 



the stars and dark matter particles within a satellite are 
central to these modern branches of astrophysics. 

In this paper, we revisit this old problem and present an 
improved analytic calculation for the tidal radius of satel- 
lites. It has been understood, at least numerically, for some 
time that the tidal radius in general depends upon four fac- 
tors: the potential of the host galaxy, the potential of the 
satellite, the orbit of the satellite and, the orbit of the star 
within the satellite. To our knowledge, however, only the 
first three factors have been addressed in anal ytic c alcula- 
tions previously (see e.g. Ivon Hoernerlll957l and lKina ll962h 
Here, we present an analytic calculation which includes the 
effect of the orbit of the star too. 

That the orbit of the sta r is important ha s been 
noted several times in the past. iToomre fc Toomrel (I1972T) 
found numerically that prograde orbi ts are more easily 
stripp ed than retrograde orbits, while iKeenan fc Innanenl 
(1975) showed that both radial and prograde orbits are 
more easily strippe d tha n retr ograde orbits. More re cently, 
iKazantzidis et all J2004) and iKravtsov et alJ ((2004) have 
found that the satellites in their simulations show tangen- 
tial velocity anisotropy near the tidal radius and that radi- 
ally biased velocity distributions are more rapidly stripped. 
Similar results have been observed in globular cluster and 
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Figure 1. A schematic diagram of the analytic set-up. The satel- 
lite and the star within the satellite both orbit in the x-y plane 
with angular velocities Q and Q a respectively. We consider the 
problem from a frame centred on the host galaxy rotating with 
angular velocity f2(t). In this frame the centre of mass of the satel- 
lite is stationary. In this diagram, the star (marked by the black 
star) is on a prograde orbit; see text for further details. 



star cluster simulation s jSeitzerjl985yTakah^ 

iGiersz fc Heggiell997l and lBaumeardt fc Makinol200.4 and 

in numerical studies of the restricted three-body problem 
jHenorJll970l) . All of these numerical results point to the 
fact that the orbits of stars play an important role in deter- 
mining the tidal radius for the satellite. In this paper, we 
shed some analytic insight on these numerical observations. 

This paper is organised as follows: in section we 
present our calculation for the tidal radius of spherical sys- 
tems. We include the effect of the orbit of the star within 
the satellite by considering three limiting cases of interest: 
prograde and retrograde circular star orbits, and pure ra- 
dial star orbits. In section [3] we then compare our analytic 
formulae with detailed N-body simulations of tidal stripping 
and show that they provide an excellent fit to the simulation 
data. Finally, in sectional we present our conclusions. 



2 AN ANALYTIC CALCULATION OF THE 
TIDAL RADIUS 

A schematic diagram of the analytic set-up is shown in Fig- 
ure Q We analyse the problem from a frame rotating with 
angular velocity, Q(t), centred on the host galaxy. We can 
choose Q(t) such that the centre of mass of the satellite will 
always b e stationary in this frame. Using Lagrangian meth- 
ods as in lBinnev fc Tremaind dl987t) . the equation of motion 
for the centre of mass of the satellite in this frame is given 
by: 



(1) 



where x is the vector distance joining the centre of mass of 
the satellite to that of the host galaxy, and V$ 9 (x) is the 
force per unit mass due to the host galaxy potential, & g (x), 
on the centre of mass of the satellite. 

In the same frame as above, the equation of motion for 
a star at an instantaneous distance, x s , from the host galaxy 
is given by: 



where V$ a (r t ) is the force on the star due to the satellite 
and r t — x s — x is the instantaneous tidal radius of the star. 
Note that from here on we will assume in this calculation 
that the star is at its instantaneous tidal radius, r = r t . 

For the special case of stars on pure circular or pure ra- 
dial orbits (two limiting extremes of interest), we may write: 



= x + n s A r t 



(3) 



where Q g is the angular velocity of the star about the centre 
of mass of the satellite. Equation is valid for circular 
orbits because in this case the space velocity of the star is 
due only to the velocity of the satellite galaxy plus that 
due to its rotation, £2 s . It is also valid for pure radial orbits 
because in this case, not only does Q s — > 0, but crucially, 
at the apocentre of the star's orbit, the radial component 
of the star's space velocity about the satellite must also be 
zero. Since all stars on radial orbits with apocentres larger 
than the instantaneous tidal radius will be stripped, we may 
take the limiting case where radial star orbits at the tidal 
radius are at apocentre. 

Equating equations Q and ©, substituting for equa- 
tion and using the fact that at the tidal radius, r t = 
x_ 3 -— x = 0, we obtain the following equation for rt: 

F. + F = 



-/ 

where 

F_ = V$ s (x) — V$ 9 (x s ) - V$ s (r t ) 

F f = A r t - A (Q A r t ) - 2f2 A (fi s A r t ) 



(4) 

(5) 
(6) 



Notice that the Coriolis terms, 2fiAi, have cancelled. For the 
case of stars on pure radial orbits where fi 3 — > 0, all Coriolis 
terms disap pear which is why in the original calculation by 
iKind <ll962T) . such terms can be left out. 

Equations (jj and © are still quite unwieldy and give 
a continuum of tidal radii for stars on different orbits and 
different alignments to the orbital plane of the satellite. We 
may greatly simplify things, however, by searching for the 
smallest tidal radii. This occurs for the co-planar orienta- 
tion shown in Figure Q By placing the star in the same 
orbital plane as the satellite, the star's rotational velocity 
then maximally adds to (prograde orbits) or subtracts from 
(retrograde orbits) the space velocity of the satellite about 
the host galaxy. Furthermore, we are interested in the phase 
of the star's orbit which similarly minimises the tidal radius. 
This is when the star is aligned along the line of centres of 
the host galaxy and satellite as shown in Figure At this 
point the star is both at its closest approach to the host 
galaxy and at the point where its rotational velocity is of 
maximal effect. 

Using the geometry of Figure^and taking just the force 
component along the line of centres of the host galaxy and 
satellite (along x), we obtain: 



F f ■ x + F ■ x = 



where 

F-x=^ 
dx 



dx 



\x — r t \ 



d$ s 
dr 



• x — — Q, r t — 2aD,Q, s r t 



(7) 

(8) 
(9) 
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and 



1 prograde 

radial 

— 1 retrograde 



(10) 



Notice that the terms involving Q_ vanish in this geometry 
and that we have introduced the notation: f2 = and 
similarly for the other vectors. 

For the extremal cases of pure circular and pure radial 
orbits, we now have three tidal radii for prograde, radial and 
retrograde orbits. It is instructive to rearrange equation Q 
to give: 



GM g (x) GM g {x-n) GM s (r t ) 



+- 



-n 2 r t -2aQ,n s r t = 0(11) 



where M a (r) and M g (x) are the satellite and host galaxy 
mass distributions respectively. Two important points can 
be made from equation 11 U . First, the fictitious forces will 
be maximised for retrograde orbits and minimised for pro- 
grade orbits. Thus, we expect that prograde orbits are more 
easily stripped than radial orbits; while radial orbits are 
more easily stripped than retrograde orbits. This will be 
tested against numerical experiments in section |S| Secondly, 
if there were only one tidal radius (as is the case for Q s =0), 
then there will always be point mass solutions for M g (x) 
and M s (r t ) which recover this tidal radius. This is the fa- 
miliar result that a single tidal radius contains information 
only about the total enclosed mass of the satellite and host 
galaxy. This is no longer true once the star orbits are taken 
into account. In this case Q, s (rt) 2 = GM s {rt)/r'i 7^ and 
the three limiting tidal radii will depend on the mass distri- 
bution of the satellite as well as the total enclosed mass. This 
point is a direct and important consequence of including the 
effect of star orbits in the calculation of the tidal radius. 

For spherical symmetry, the energy equation for the 
satellite may be rearranged to give: 



J 2 = 



E = 



2x1 



(12) 



(13) 



where J is the specific angular momentum of the satellite, 
E is the specific energy and x p and X a eLT6 the satellite peri- 
centre and apocentre respectively. The satellite orbit is then 
fully specified by x p < x < x a and we may readily obtain 
fl(t) = J/x(t) 2 from equation ltT2t . 

Since £1 = Q(t) and x = x(t), for all but circular satel- 
lite orbits, the three limiting tidal radii solutions to equation 
1111 are a function of time: r t = r t (t). Whether a time in- 
dependent value for the tidal radius is appropriate in such 
situations is discussed further in section EP1 

2.1 Point mass potentials 

Given <& g (x), $ s (rt), x(t), x v and x a , we now have enough 
information to solve equation jilt for the three values of 
rt(t). In general, this solution must be found numerically. 
An analytic solution may be obtained, however, for the spe- 
cial case of point mass potentials for the satellite and host 
galaxy, in which case the satellite's orbit is Keplerian. In 
this case, equation 1121 gives us: 



9/ 



GM a 



a 2x a Xp , 2\ 

A = — = a(l — e ) 



(14) 
(15) 



where G is the gravitational constant, M g is the mass of the 
host galaxy and we define here the standard Keplerian semi 
major axis, a, and eccentricity, e, as: 

e = Xa ~ x P (16) 



a ~\~ Xp 



(17) 



In fact, these definitions involve only x a and x p , and so hold 
good for any orbit in a spherical potential. 

Substituting equation I14H into equation J7J, using 
point mass potentials, and solving, we obtain: 



r t 



x^A 



M a 



1/3 



^JcP- + 1 + 2x/A - 
A + 2x 



2/3 



(18) 



galaxy respectively and a is as in equation 1101 

Equation (|180 gives rt(x) or si nce x = x(t), equival ently, 
rt(t). As o riginally suggested by Ivon HoerneJ (119571) and 



King ( 1962), we may take the limiting case which gives the 
smallest tidal radius - namely that at pericentre; x — x p . In 
this case, equation (1181 reduces to: 



r t 



Ms 



1/3 



1 



1/3 



la 2 + 1 + 



2/3 



1 



2 

l + e 



(19) 



where e is the eccentricity of the satellite orbit. 

It is reassuring to note that, for the case of pure radial 
star orbits, a — and equation 119H reduces to the familiar 
tidal radius: 



rt 



M s 



M 9 (3 + e) 



1/3 



(20) 



The above analytic forms for point masses are useful 
in explicitly showing that r t has three limiting values at 
each position, x p < x < x a , of the satellite. From equation 
I19H . we see that the prograde stripping radius (a = 1) is 
smaller than the radial one (a — 0) which is smaller than 
the retrograde one (a — —1), as expected. 

2.2 Power law potentials 

In practice, point mass potentials are a poor approximation 
to observed galaxy mass distributions which are much more 
extended. Another set of analytic solutions to equation 1111 
may be obtained for power-law density profiles of the form 
ps t3 — A Stg r~"' s - a , where as previously s and g subscripts de- 
note the satellite and galaxy respectively 1 . Solving equation 
II II for these density profiles gives: 



4 s A'(3- 7 a) 



4i 4 



a + 



7s a; 



3- 79 A' 



+ 1 



-2/7= 



(21) 



1 Such a den sity profile is only p hysical for a restricted range of 

2 ^ 7s, 9 < 3 iBinnev fc Tremainelll987l) . 
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where 



A' = 2 



X a )x n Xf, 



(3 - 79 )(2 - 79 )(4-x2) 



(22) 



In the limit x p — x a (a circular satellite orbit), A' = ^z-i 
and so is well-behaved. 

A particular case of interest is for isothermal sphere 



density profiles where 7 S 



2. This is because these ap- 



proximate well the observed mass distributions in galaxies. 
In this special case, equation 1211 reduces to the following 
simple form: 



n 



y/k"A s /A g -a+ ^a 2 + l + x 2 /A" 



■A"/x 2 



where 

A" = 2 f aX \ In ^ 

>*jq J^a 

In the limit Xv - 



(23) 



(24) 



A" = x%. 



For the isothermal sphere, A s /A g is more usually writ- 
ten a 2 /a g , where a s>g is the velocity dispersion of the 
satellite/host galaxy. For satellites on pure circular orbits 
(x a — x p ), with stars on pure radial orbits (a = 0), equa- 



tion!23lrcduces to: rt 



x/2 CT „ 



, in agreement with the formula 



derived recently bv lDiemand et all i2005l) . 

Another interesting limit to equation 12 II is for 7 Sj9 = 3. 
This limit recovers the result for point mass potentials given 
in equation 1181 



2.3 Split power law potentials 

Power law potentials provide a useful class of analytic so- 
lutions. However, all valid power law solutions except the 
point mass have diverging mass at large radii and a limited 
range of applicability. (Recall that 2 ^ 7 Sj9 ^ 3 in these 
potentials.) Split powe r law density profiles do n ot suffer 
from the se restrictio ns (lHernquisdll990t ISahalll992l iDehnenl 
ll993l and lZhaal l996). A fully analytic sub-set of general split 
power law profiles is given by: 



Ps,g 



A/ 3 , 9 (3- 7s , 9 ) 



4vrr|, s 



{r/r s , g )^(l + r/r s , g ) 4 



M 3 , g G 



- ls,g, 



[a 



,/ry 



1] 



(25) 



(26) 



, g are the mass and scale length in both 
cases, and 7 S , 9 is the central log-slope of the split power law 
density profile. 

Equation II II may be solved using the above split power 
law profiles, but is only fully analytic for 7s = and r t <C 
x + r g . The second restriction is nearly always an excellent 
approximation since it is likely that both r t < i and r± <C 
r g . The first restriction works for cored satellite halos, and 
for satellite halos with rt ^> r s for which only the outer 
regions are being sampled. In these limits, eouation llll gives: 



n 



x^A" 



1/3 



yV + 1 - qx A /A" 



2/3 



A" 



qx* 



-r s (27) 



(x + r 9 ) 3 "Ts 



1-7b + 



la 



A" 



r g (1-1g){x1,-xl) 



l + r g /x 
[(l + rg/xp 



(28) 



{l+Tg/Xa 



( \ 3 — 79 

In the limit x p = x a , A'" — x p ( — 5 — ) • I n the limit 

\ x p-r r g j 

Tt 3> r s and i»r s , equation 1271 reduces to eouation ll8l for 
point mass potentials, as expected. 

In the following section, where we compare the theoret- 
ical tidal radii with N-body simulations, we solve equation 
Ijllll numerically for the same potentials as those used in the 
simulations. Using point mass potentials instead can lead to 
errors in the tidal radii as large as a factor of two; using 
isothermal potentials would give smaller errors, since these 
well-approximate the galactic potentials used in the simula- 
tions. 



3 COMPARISON WITH NUMERICAL 
SIMULATIONS 

We set up the initial conditions for the satellite by drawing 
the positions from an analytic density profile and the ve- 
locities from a numerically calculated distribution function. 
This proc ess is described an d tested in detail in a compan- 
ion paper jRead et alJl2005l) . We note here that using distri- 
bution functions rat her than the mor e familiar Maxwellian 
approximation (as in iHernqui st 1993) for calculating parti- 
cle velocit ies is much more accura te for simulations of tidal 
stripping jKazantzidis et aT]l2004f) . 

We used a Plummer density profile for the satel- 
lite which is a self-gravitating, single component, spherical 
gal axy comprising only stars. The Plummer profile is given 
by jBinnev fc TremaindlT987h : 

_ 3M S 1 r , m 
Pplum ~ 47ra| (1 + r|)5/2 (M) 

where M s = 10 7 Mq and a B — 0.23 kpc are the mass and 
scale length of the satellite respectively. The satellite veloc- 
ity distribution was initially isotropic and when evolved in 
isolation was found to be extremely stable over a Hubble 
time. We used a simulation resolution of 10° particles for 
the satellite and a force softening of 10 parsecs; the simula- 
tions were found to be very well resolved - higher resolution 
test runs with 10 6 particles and parsec scale force softening 
produced converged results. 

The initial conditions were evolve d using a modified ver - 
sion of the GADGET N-body code JSpringel et al.l 1200 ll) 
modified to permit a fixed potential to model the host 
galaxy. We used a host galaxy pot ential chosen to pro- 
vide a good fit to the Milky Way jLaw et alJl2005l) . We 
used a Miyamoto-Nagai potential for the Milky Way disc 
and bulge and a logarithmic potential for the Milky Way 
dark matter halo. These are given by respectively (see e.g. 
iBinnev fc Tremaindll987l) : 

<3> mn (R,z)= ~ GMd (31) 



y/R 2 + (a + Vz 2 + b 2 ) 2 



where 



where M — 5 x 1O 1O M0 is the disc mass, a = 4 kpc is the 
disc scale length and b — 0.5 kpc is the disc scale height, 
and: 



(29) 
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Figure 2. The tidal stripping of satellite stars as a function of the stellar orbits. Each panel shows a different time output as the satellite 
orbits around the host galaxy on a circular orbit at 80 kpc. The blue, green and red lines show the radial distribution of stars which 
were initially on prograde, radial and retrograde orbits respectively. The solid lines show the initial distributions, while the dotted lines 
show the evolved profiles. The vertical dashed lines show the three analytic tidal radii calculated from equation [7] prograde (blue), radial 
(green) and retrograde (red). 



^iog(f) = m (Rc + r 2 ) + constant (32) 

where R c — 4.1 kpc is the halo scale length and vo = 
220km/s is the asymptotic value of the circular speed of 
test particles at large radii in the halo. 

We held the satellite properties and the host galaxy 
properties fixed and varied only the orbit of the satellite 
about the host galaxy. We use a circular orbit at 80 kpc 
(section 13. II and an eccentric orbit with apocentre 85 kpc 
and pericentre 23 kpc (section 13.31 . These orbits are typ ical 
of Local Group dS ph galaxies (see e.g. lPiatek et all2002l and 
iPiatek et ai1l2005l) . 

Finally, it is not trivial to mass and momentum centre 



the satellite when performing analysis of the numerical data. 
An incorrect mas s centre can lead to spurious density and 
velocity features jPontzen et al. 2005). We use the method 
of shrinking sp heres to find the m ass and momentum centre 
of the satellite JPower et alJl2003ft . 



3.1 Circular satellite orbits 

The first numerical test is a circular orbit of a satellite which 
experiences relatively weak tidal effects from the host galaxy. 
This is because the analytic formulae assume that the satel- 
lite potential remains constant. If a large fraction of the 
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satellite's mass is stripped away, this approximation starts 
to break down. In this case, the tidal radii should be recalcu- 
lated at each time step iteratively as mass is removed, much 
as i s done in semi-analy tic simulations of tidal stripping (see 
e.g. iTaffoni et ai]l2003T) . 

We wish to understand how a satellite star is affected 
by the host galaxy's tidal field, as a function of the star's 
orbit. In order to do this, we calculate the orbit of each star 
in the satellite initially. For a spherical system, the energy 
equat ion for an individual star can be simply rearranged to 
give feinnev fc Tremainelll987l) : 

r 2 = 2{E s -<f> s (r))- J s 2 /r 2 (33) 

where r is the radial velocity of the star with respect to the 
satellite's centre of mass, E s is the star specific energy, J s is 
the star specific total angular momentum and $ s (r) is the 
potential of the satellite. 

At pericentre, r p , and apocentre, r a , r = and so we 
may calculate r a and r p for each satellite star by finding the 
roots of equation (I33H . We may then obtain the eccentricity 
of each star orbit, e, as defined in equation 1161 . In order 
to do this, we must know $ s (f)i and J s and E s for each 
star. The specific kinetic energy and angular momentum for 
each star may be simply computed from their phase space 
coordinates. For 10 star particles, the potential due to the 
self gravitating stellar distribution is prohibitively slow to 
compute via direct summation. To avoid this problem, we 
use a tree code to compute the potential foehnerJl2000l) . 

Ideally, we would like to compare pure circular and pure 
radial orbits with the analytic calculation in section [5] How- 
ever, models built either with pure circular or pure r adial or- 
bits are unstable (see e.g. lBinnev fc Tremainell98'it) . We use 
instead an isotropic distribution of velocities which has the 
advantage that it is stable, but the disadvantage that very 
few stars are on either pure circular or pure radial orbits. 
To obtain good statistics, we define an orbit as being 'cir- 
cular' if e < 0.5 and 'radial' if e > 0.7. The 'circular' orbits 
are then divided into prograde ( J s • J > 0) and retrograde 
( J_ s • J < 0) orbits; where Js_ and J are the star and satellite 
specific angular momentum vectors respectively. The radial 
orbit boundary (e > 0.7) is chosen such that there are as 
many stars on 'radial' orbits initially as there are on pro- 
grade or retrograde orbits. 

Having divided up the initial stellar distribution in this 
way, we can track the evolution of prograde, radial and ret- 
rograde stars as the satellite orbits around the host galaxy. 
Figure shows four snapshots in time over 10 Gyrs as the 
satellite orbits around the host galaxy. 

From Figure^ we can see that the analytic model works 
very well for the first Gyr. The dotted lines show that strip- 
ping occurs for each of the orbit-types only up to the an- 
alytic tidal radii. Note that we have included both bound 
and unbound stars in this plot. This is why the dotted lines 
turn upwards at large radii. It is a well-known result that 
unbound stars can take a long time to actually drift away. 
However, that is a different effect to the one being discussed 
here. The retrograde stars interior to the retrograde strip- 
ping radius are bound energetically and not just unbound 
stars taking a long time to escape. We have explicitly calcu- 
lated energies to check this. 



Over longer time scales 2 , the analytic limits appear to 
fare less well: all of the orbit types begin to become stripped 
into the prograde limit (blue dashed line). The natural ques- 
tion is then, what is causing this effect? Two obvious pos- 
sibilities are that either the potential of the satellite has 
changed significantly, in which case the tidal radii should 
migrate inwards, or that the orbits of the stars have been 
altered by the tidal field of the host galaxy. By plotting the 
potential of the satellite as a function of time, it is easy to 
show that option one is not correct. As might be expected 
given the small number of stars which are actually stripped 
in this simulation, the satellite potential remains almost con- 
stant over the whole simulation time of 10 Gyrs. Rather, the 
orbits of the stars are changing in the tidal field of the host 
galaxy. We discuss this further in the following section. 

3.2 Star orbits in a tidal field 

So far we have implicitly assumed that the orbits of stars 
within the satellite will be unchanged by the tidal field. This 
is not the case, and in fact determining the nature of such 
orbits within a tidal field is an interesting problem in its 
own right which is beyond the scope of this present work. 
However, what matters for determining the tidal radius is 
only the relative velocity of the star with respect to the cen- 
tre of mass of the satellite. In this sense, our three limiting 
tidal radii for prograde circular, pure radial and retrograde 
circular orbits delimit regions of phase space which, if a star 
ever enters, then it will become stripped. 

In Figure |3 we plot two star orbits taken from the sim- 
ulations. If the star orbits were evolved in isolation from a 
tidal field, then the star angular momentum would be con- 
served (since the satellite potential is spherical). In a tidal 
field, the star angular momentum, J s , is not conserved. The 
z-component, J s ,z, can alter, causing a flip from retrograde 
to prograde motion. Such stars can thus become stripped 
even if they start out on orbits which, if evolved in tidal iso- 
lation, would not be. Stars most likely to be affected in this 
way are those with low J S;Z - those on near-polar, or highly 
eccentric orbits. 

In the first orbit (top three space projections in Fig- 
ure tidal forces act along the direction of motion of the 
star. In this case, an initially prograde orbit loses its tangen- 
tial velocity and is converted to a radial orbit. After further 
deceleration, the final orbit is only very slightly retrograde 
and appears nearly radial in the plot. In the second orbit 
(bottom three space projections), the star is on a more po- 
lar orbit. Tidal forces acting perpendicular to star's orbit 
plane tilt the orbit about the z-axis, causing a flip from ret- 
rograde to prograde motion. Both orbits have low J s , z : the 
first because it is eccentric; the second because it is polar. 

The result of these orbital transformations is that the 
three tidal radii of section [2J over long times, converge on 
the prograde stripping radius. This is because any star be- 
yond the prograde stripping radius which has its orbit trans- 
formed to be prograde is rapidly stripped. Over long times 
it becomes increasingly likely that low angular momentum, 

2 The time for one circular orbit at 80kpc in the potential we 
use is ~ 2 Gyrs. Throughout this paper we refer to 'short' times 
as being ,$ 2 Gyrs and 'longer times' as being 2 Gyrs. 
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high eccentricity, stars will for a short time be pushed onto 
prograde orbits, at which point they will be carried away by 
tides. 

3.3 General satellite orbits 

In this section, we consider the more complicated case of a 
satellite on a highly eccentric orbit. We use an orbit with 
apocentre 85kpc and pericentre 23kpc. FigureQ]shows four 
snapshots in time over 10 Gyrs. 

As discussed in section |5J the three limiting tidal radii 
for prograde (blue dashed line) , retrograde (red dashed line) 
and radial (green dashed line) are functions of time. This 
is why the solid lines are in different places in the panels 
of Figure 2] The dashed vertical lines in Figure 2] show the 
values of the three tidal radii at the pericentric distance of 
the satellite; these do not vary with time. An added compli- 
cation is that the galactic potential we use is not spherical. 
As a result, the orbit of the satellite will precess and its 
angular momentum vector, J, will change direction. Thus 
defining 'prograde' and 'retrograde' stars becomes formally 
more difficult as such terms now become a function of time 
along the satellite's orbit. With the simulation parameters 
chosen here, however, the precession of the satellite's orbit 
plane is small over 10 Gyrs and makes little difference to the 
initial orbit classification. Because of this, we use the initial 



value of J when defining prograde, radial and retrograde 
stars. 

Notice that there are now two points at which the satel- 
lite stars show a tidal 'break radius'. One is at the current 
tidal radius, the other is close to the pericentric tidal ra- 
dius. Similar results for the surface brightn ess distributions 
have been observed by other authors (see e.g. l.Tohnston et all 
2002). Each of these break radii are split into three further 
limiting tidal radii as a function of the star orbits. 

Over long timescales ( £ 8 Gyrs) , the break radius 
which lies at the current tidal radius of the satellite's or- 
bit disappears. This is because, after the initial distribution 
has been stripped of all stars beyond their tidal radii, fur- 
ther stripping can only occur as stars heated by the tidal 
field of the host galaxy migrate outwards. This effect is much 
smaller and produces only a tiny perturbation on the surface 
brightness profile which, at these late times, now appears 
much smoother. 

Notice that over a Hubble time, the central density of 
all orbits lowers in its normalisation, even in the very centre 
of the satellite. This is not due to tidal stripping, but rather 
to a different physical effect which occurs only for highly 
eccentric orbits, or for orbits which cause the satellite to 
pass thr ough the plane of the disc - namely tidal sho cking 
fsee e.g. lOnedin OstrikerllTOol. lOnedin et aUliflOfli and 
iGnedin et alJll999l) . These tidal shocks can (and in this case 
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Figure 4. The tidal stripping of satellite stars as a function of the stellar orbits. Lines and panels are as in figure|2] but now the satellite 
is on a highly eccentric orbit with apocentre 85 kpc and pericentre 23 kpc. The vertical solid lines show the three tidal radii calculated 
at the current time; while the vertical dashed lines show the three tidal radii calculated at the pericentre (distance of closest approach) 
of the satellite's orbit. 



do) become more important than stripping in the very centre 
of the satellite. We discuss the effects of tidal sh ocks in more 
detail in a companion paper, jRead et aUl2005ft . 



As for the circular orbit case, orbital transformations 
drive all star orbits eventually towards the prograde limit at 
pericentre (blue solid vertical line) . However it is difficult to 
ascertain whether the prograde limit is ever actually reached 
in practice for general satellite orbits because of the action 
of tidal shocks. 



3.4 Implications of the new tidal radii 

A natural question relates to the observability of the three 
tidal radii in the projected surface brightness profile and 
velocity dispersion of the stars. These are shown in Figure 
|S] Notice that if only surface brightness information were 
available, the tidal radius determined purely from 'features' 
in the light profile would be erroneous. For the case of a 
satellite on a circular orbit (top panels), the surface bright- 
ness profile is depleted at the prograde stripping radius (blue 
vertical line) and breaks at the retrograde stripping radius 
(red vertical line) . For a more general satellite orbit (bottom 
panels) there are no observable features in the light profile at 
any of the analytic tidal radii. This is because tidal shocking 
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Figure 5. The projected surface brightness distribution (left) and velocity dispersion (right) for a satellite on a circular orbit (top panels 
and see section 13.11 and an eccentric orbit (bottom panels and see section 13.31 . The black lines show the initial profiles while the red 
lines show the profiles after ~ 10 Gyrs. The solid vertical red, green and blue lines show the three analytic tidal radial from equation 171 
In the plot of the velocity dispersions, the solid red line shows the radial velocity dispersion while the dotted and dashed lines show the 
9 and 4> components. 



washes out such information (notice that such tidal shocks 
also lower the central surface brightness as the satellite be- 
comes puffed up). However, if the radial and tangential ve- 
locity dispersions are known then the tidal radii may be 
much better determined. Notice that for both the satellite 
on the circular orbit and the more general orbit, tangential 
velocity anisotropy appears at the prograde stripping radius. 
This is because, as discussed in section 13.21 stars on radial 
or near radial orbits are, over long times, pushed momen- 
tarily onto prograde orbits, at which point they are stripped 
away. Thus even at the prograde stripping radius, there is a 
depletion of radial orbits with respect to retrograde orbits 
and this leads to the onset of tangential anisotropy. Such 



anisotropies increase outwards as stars moving on radial or- 
bits become more easily stripped at the radial stripping ra- 
dius. The presence of tangential anisotropy provides a much 
more robust signature of the true tidal radius than can be 
measured from the light profile alone. 



4 CONCLUSIONS 

We have presented an improved analytic calculation for the 
tidal radius of satellites and tested our results against N- 
body simulations. 

The tidal radius in general depends upon four factors: 
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the potential of the host galaxy, the potential of the satel- 
lite, the orbit of the satellite and the orbit of the star within 
the satellite. We demonstrated that this last point is criti- 
cal and suggest using three tidal radii to cover the extrema 
range of orbits of stars within the satellite. In this way we 
showed explicitly that prograde star orbits will be more eas- 
ily stripped than radial orbits; while radial orbits are more 
easily stripped than retrograde ones. This result has previ- 
ously been established by several authors numerically. We 
showed further that, in general, the tidal radius does not 
depend only on the enclosed mean density of the satellite. 
While for stars on pure radial orbits it is reasonable to ap- 
proximate the tidal radius as being that where the density 
of the satellite matches that of the host galaxy, for prograde 
and retrograde star orbits, the tidal radius also depends on 
the mass distribution of the satellite galaxy. 

Over short times ( ,$ 1 — 2 Gyrs ~1 satellite orbit), we 
find excellent agreement between our analytic and numerical 
models. Over longer times, star orbits within the satellite 
are transformed by the tidal field of the host galaxy. In a 
Hubble time, this causes a convergence of the three limiting 
tidal radii towards the prograde stripping radius. 

Satellites which have been tidally stripped will show 
tangential velocity anisotropy near their tidal radii, with a 
depletion of both prograde and radial orbits relative to the 
retrograde orbits. This must be true irrespective of the ini- 
tial conditions. These results na turally explain the numer - 
ical observations dating back to iToomre fc Toomrd (ll9T2f) 
and iKeenan fc Innanenl jl975T) which found exactly these 
velocity anisotropies in their simulations of galaxy-satellite 
interactions. 

In the future, with improved kinematic data for nearby 
satellites from SIM, it may be possible to search for tan- 
gential velocity anisotropies as a 'smoking gun' from tidal 
stripping, and thereby much better determine the tidal radii 
of nearby galaxies and star clusters. This would represent 
a significant improvement on current methods which have 
historically found 'tidal radii' from features in the surface 
brightness profiles alone. 
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